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ABSTRACT 

We present ground-based SpectroCam-10 mid-infrared, MMT optical, and 
Spitzer Space Telescope IRS mid-infrared spectra taken 7.62, 18.75, and 19.38 
years respectively after the outburst of the old classical nova QU Vulpeculae 
(Nova Vul 1984^2). The spectra of the ejecta are dominated by forbidden line 
emission from neon and oxygen. Our analysis shows that neon was, at the first 
and last epochs respectively, more than 76 and 168 times overabundant by num- 
ber with respect to hydrogen compared to the solar value. These high lower limits 
to the neon abundance confirm that QU Vul involved a thermonuclear runaway 
on an ONeMg white dwarf and approach the yields predicted by models of the 
nucleosynthesis in such events. 
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Subject headings: circumstellar matter: classical novae, cataclysmic variables - 
stars: individual (QU Vul, Nova Vul 1984#2) 



Introduction 



Gehrz et al.l (119851 ) discovered a very strong [Ne II] 12.8 forbidden emission line in the 
classical nova QU Vulpeculae (Nova Vul 1984#2) 140 days after its outburst. The presence 
of this strong neon line and the subsequent appearance of additional infrared metallic forbid- 



den e mission lines in the spectrum as QU Vul evolved (IGehrz et al.lll986l ; [Greenhouse et al. 



19881 ) confirmed that there are oxygen-neon-magnesium (ONeMg) white dwarfs (WDs) in 
classical r iova binary systems whose erup t ions could significantly e nr ich the local i nterst ellar 
medium jGehrzl Il988l : iGehrz et al. l ll998l : lGehrzlll999l . bood 120081 ) . iGehrz et all Jl998h ar- 
gued that theoretical simulations of nova thermonuclear runaways (TNRs) provide evidence 
that ONeMg novae may, in fact, compete with supernovae as sources of some CNONeMgAl 
isotopes. We have pursued a long term program to measure metal abundances in nova ejecta 
using infrared (IR) spectroscopy of forbidden emission lines to test this hypothesis. 

Measurements of the integrated fluxes of these electron-impact (coUisionally) excited, 
forbidden lines in nova ejecta many years after outburst are particularly valuable for eval- 
uating abundances. In this case, the density of the expanding ejecta will have fallen well 
below the level where the lines are quenched by collisions, and the line fluxes are entirely due 
to spontaneous emission. Since coefficients of spontaneous emission are fairly well known for 
many important IR lines, the absolute numbers of atoms involved in these transitions can 
be calculated with reasonable certainty given a known distance. Because the ejected mass of 
hydrogen is well known from early infrared and ra dio emission measures when most of the 
gas is ejected and when it is almost entirely ionized (lGehrzlll999l . |2002| ) . the ratio of the mass 
of an element required to produce the observed IR forbidden emission fines to the ejected 
hydrogen mass gives a fairly well constrained lower limit to the abundance of that element. 
When emission from all the lines from all ionic species and all isotopes of an element can be 
observed and evaluated, the lower limit converges to the actual abundance. 

Several mid-IR neon forbidden emission lines in QU Vul are sufficiently strong that we 
have been able to observe them using SpectroCam-10 on the 5 m Hale Telescope at Mount 
Palomar and the Infrared Spectrograph (IRS) on the Spitzer Space Telescope at epochs of 
7.62 and 19.38 years after the eruption. These measurements have enabled us to determine 
a fairly strong lower limit to the neon abundance in this nova as reported below. 
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OBSERVATIONS AND REDUCTION 



2.1. Ground-based Mid-IR 

We observed QU Vul with the 200 inch Hale Telescope at Mount Palomar0 on 1992 Au- 
gust 09.3 UT (day 2783.9, 7.62 years post-erup tion) using the SpectroCam-10 (SC-10) mid- 
IR camera/spectrograph (IHayward et al.lll993l ) in its long-slit, low-resolution spectroscopy 
mode with a 2" wide slit and a resolving power of R = A/AA ~ 150 at 12.81 /im. We note 
that the entire extent of the ejected shell was contained within the slit (see Krautter et al. 
2002). A single grating setting provided a 10.6 to 13.2 fim spectrum that included the [Ne 
II] 12.81 /im line. The total on-source integration times were 160 s on QU Vul and 20 s on 
the standard star a Lyr. The data were reduced by extracting 1-dimensional spectra from 
the 2-dimensional spectral images, dividing the n ova by the st andard, then multiplying by 
the absolute spectral energy distribution of a Lyr ( 1Cohenlll992l ). The calibrated result, with 
2-pixel binning for clarity, is shown in Figure [H In Table [T|, we give the total flux measured 
in the line. 



2.2. Ground-based Optical 

An optical spectrum of QU V ul was obtained with the 6.5 m Multiple Mirror Telescope 



and Blue Channel Spectrograph (ISchmidt et al.lll989l ) on 2003 September 13.31 UT (day 



6836; 18.72 years post-eruption) in generally photometric conditions and 1" FWHM seeing. 
Three 600 s integrations of the nova were obtained through a 1" x 180" entrance slit oriented 
along the parallatic angle to avoid loss of light at the slit from differential atmospheric re- 
fraction and to improve the relative photometric accuracy. A 300 lines/mm grating centered 
at 6400 A was utilized which covered the entire optical region on the detector; however, a 
long-pass order separation filter blocked all light shortward of ~3800 A and yielded cali- 
brated spectra covering the region 3800 - 8000 A at a resolution of about 7 A FWHM. The 
spectrophotometric standard star BD+28 4211 was observed at a similar airmass as QU Vul 
using a wider 5" x 180" entrance slit to enable accurate flux calibration. The spectra of HeN- 
eAr and a quartz-halogen lamps were obtained to facilitate accurate wavelength calibration 
and fiatfielding respectively. The spectra were reduced using standard IRAF package^ and 



^Observations at the Palomar Observatory were made as part of a continuing collaborative agreement 
between the California Institute of Technology and Cornell University. 

^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
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spectral extraction techniques. The optical spectrum is presented in Fig. [2l 



2.3. Spitzer IRS 

We obse rved QU Vul on 2 004 May 11.7 UT (day 7077.6, 19.38 y ears post-eru 



ption) 



2004 



with the IRS (IHouck et al.ll2004f) on the Spitzer Space Telescope (Spitzer. \Werner et al. 
Gehrz et al.l 120071 ) as part of the Gehrz Guaranteed Time Observing Program (GGTOP), 
program identification number (PID) 124. The observations were conducted using the short 
wavelength (5.2 - 14.5 /xm) low resolution module (SL), the short wavelength (9.9 - 19.6 /xm) 
high resolution module (SH), and the long wavelength (18.7 - 37.2 fim) high resolution 
module (SH). The resolving power for SL is i? = A/AA ~ 60 — 120, while R ~ 600 for SH and 
LH. Emission lines observed with the latter modules are marginally resolved. All observations 
utilized visual PCRS peakup on nearby isolated stars to ensure proper placement of the 
target in the narrow IRS slits. The spectroscopic astronomical observation request (AOR) 
key 0005044992, consisted of 5 cycles of 6 second ramps in SL (30 seconds on-source), and 6 
cycles of 6 second ramps in SH and LH (36 seconds on-source each). 

IRS Basic Calibrated Data (BCD) products were calibrated and processed with the 
Spitzer Science Center (SSC) IRS pipeline version 13.2. Details of the calibration and raw 
data processing are specified in the IRS Pipeline Description Document, vl.olfl Bad pixels 
were interpolated over in individual BCDs using bad pixel masks provided by the SSC. 
Multiple data collection events were obtained at two different positions on the slit using 
Spitzer 's nod functionality. Sky subtraction was only possible for the SL observations as 
no dedicated sky observations were performed for the SH and the LH mode observations. 
Sky subtraction was performed by differencing the two dimensional SL BCDs to remove the 
background flux contribution. Spectra were then extracted from the background corrected 
SL data and the SH and LH BCDs with SPICE (version 1.3-betal) using the default point 
source extraction widths. The extracted spectra were then combined using a weighted linear 
mean into a single output data flle. The high resolution IRS data were not defringed. At the 
time of reduction, the errors generated by the SSC pipeline were not reliable and therefore 
errors were estimated from the standard deviation of the flux at each wavelength bin. The 
continuum was not detected over much of the spectral range of our data (absolute flux 
calibration of the IRS modules is < 10%) and no emission was detected in the second order 
SL spectra. A good flt to the spectral lines was obtained using a non-linear least squares 



Science Foundation. 

•^http:/ /ssc. spitzer. caltech.edu/irs/dh/PDD.pdf 
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Gaussian routine (the Marquardt method, iBevington &: RobinsonI Il992l ) that fit the fine 
center, fine amphtude, fine width, continuum amphtude and the slope of the continuum. 
The spectra are shown in Figs. [3] through O 

Given that the [Ne II] 12.81 fim line is accessible from the ground, we have investigated 
the possibility of continuing to monitor it in QU Vul using today's very large ground-based 
telescopes. During the interval between our SpectroGam-10 measurement of QU Vul (7.62 
years after outburst) and our Spitzer measurement (19.32 years after outburst), we at- 
tempted unsuccessf ully several times to det ect this line in QU Vul using the long wavelength 
spectrometer LWS (jjones fc Puetterlll993l ) on the Keck I 10 m telescope. More advanced 
mid-IR spectrometers are now available. One example is the VLT Imager and Spectrometer 
for mid-IR (VISIR) at ESO Paranal. We have used the VISIR Exposure Time Estimator 



(http://www.eso.org/observing/etc/ ) to calculate that VISIR on the ESO Paranal UTS 8.2 
m telescope would require about 15 hours of integration time to achieve the same signal-to- 
noise ratio (20/1) on this line as we reached with Spitzer IRS SL in 30 seconds. We conclude 
that further ground-based monitoring of the [Ne II] 12.81 fim line in QU Vul is untenable. 



3. DETERMINATION OF ABUNDANCES 



Lower limits to the abundance of an element in the ejecta of a nova can be determined by 
dividing the number of atoms required to produce the lines of the element that are observed 
by the total number of hydrogen atoms in the ejecta. The lower limit so determined converges 
to the actual abundance when all lines of all ions and isotopes can be evaluated. At very late 
epochs, hydrogen emission lines and the free-free continuum are no longer observed in the IR, 
so that the number of hydrogen atoms must be determined from the mass known to have been 
ejected in the eruption. We note that the majority of the mass is expelled within the first 
few hours of the outburst and very little thereafter, so that the number of hydro gen atoms is 



well constrained by early IR measurement s and radio continuum observations (jGehrzlll988 
Gehrz et al.lll998l : lGehrzlll999l . bood boosh . 



The neon and oxygen emission lines observed in the spectra of QU Vul are forbidden 
transitions excited by electron impact. The total number of ions required to produce the 
observed line strengths, assuming a uniform temperature distribution and density along the 
line of sight in a homo geneous me d ium, can be determined following detailed balancing 
arguments described by lOsterbrockl (119891 ). The [Ne II], [Ne VI], and [O IV] lines can be 
treated as arising from two-level atoms. [Ne III] has a ^P ground term giving rise to emission 
lines at 15.56 and 36.02 /xm. Our spectra do not extend to long enough wavelengths to 
measure the strength of the 36.02 fim line. In this case, we must treat the 15.56 fim line as 
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coming from a two-level atom, an approximation that gives us a lower limit to the number 
of Ne III atoms in the ejecta. 

The level populations (total number of ions), Ni and Nj, in the lower (i) and upper 
(j) energy states can be computed from the line luminosity, = AttD'^F, where D{cm) is 
the distance to the nova and F (erg s~^ cm ~^) is the integrated apparent intensity of the 
line. For QU Vul, we take D to be 3.14 kpc (jKrautter et al.ll2002l ). We regard this distance 
estimate as being superior to others given in the literature because the shell was well resolved 
by the HST NICMOS imager allowing a well constrained expansion parallax to be derived. 
The total number of ions in the state Nj is 



Lc (l + K/ricrit]) 



where Ue (cm ^) is the electron number density, Ucrit is the critical density of the transition. 



h is the Planck constant, and Uji (Hz) is the frequency of the transition. The parameter 
Qij (cm^ s~^) is the excitation rate coefficient, which is a function of temperature, and is 
expressed as 



8.63 X 10- 

UJiTe 



(2) 



where uoi is the statistical weight of the lower state, k is the Boltzmann constant, Te (K) 
is the electron temperature, and Tj .,- is the dimensionless thermally-averaged effective colli- 



sion streng th (IHummer et al. 
taken from 



1993r) . Values of T,;,, which i s temperature dependent, were 



Pradhan fc Pend (llQQsh . ISaraph k TuUvl (11994h . McLaughlin fc Belli fl200oh . and 



Blum fc PradhanI (119921 ) for [Ne VI] (^Pg/g p^^^), [Ne II] P3/2), [Ne III] 



^Pi ^3 P2), and [O IV] (2p 



3/2 



P1/2) respectively. 

The criti cal density, rirrit (cm~'^), for a given transition is computed from the ratio 
of Aij to qji (lOsterbrockl Il989l ) where Aij (s~^) is the radiative transition probabihty and 
qji (cm^ s~^) is the temperature-dependent, de-excitation rate coefficient defined as 



8.63 X 10- 



1/2 



-T 



(3) 



where ujj = (2 J + 1) is the statistical weight of the upper state and J is the total angular 
momentum quantum number of the upper state. 



When He << Ucrit, coUisional de-excitation of the ion in state (j) is not significant 
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compared with radiative decay. As can be seen from Table 1 (columns 8 and 9), this was 
the case for the ejecta of QU Vul when the Spitzer observations were made in 2004. Aji 
values obtained from the NIST databasfl are 2.02 x lO'^ s-\±10%) for [Ne VI], 8.59 x 
10-3 s-i (±10%) for [Ne II], and 5. 19 x 10-^ s-^ (±3%) for [O IV]. For [Ne III], we adopt 



5.84 X 10-3 s-i (±10%) from lKramida fc Navel J2006h . 



The total number of ions in state A^,- is expressed as 



N. 



(4) 



The electron number density, fig, of the ejecta at the epoch commensurate with the 
given neon line observation was inferred by adopting a free-expansion model for the ejected 
shell. Assuming spherical symmetry and a uniform shell with a filling factor of unity. 



3 M 



M, 



Att i?3 (Vo At)3mp 



(5) 



where Mej = 3.6 x 10-^ Mq JTavlor et al.lll988h is the mass of the e jected shell of radius 



-R.s (cm), Vn is the outflow velocity in cm s" taken to be 1190 km s'^ f lKrautter et al.ll2002l 



Rosino et al.lll992l ). At is the time since maximum light (day = 1984 December 25.1 UT 
= JD 2,446,059.6) in seconds, and rrip is the mass of the proton. The slow velocity chosen 
here is consistent with the line widths, which are marginally resolved by IRS SH and SC-10. 
We emphasize that this is the mean electron density inferred from the unity filling factor 
and that there is some evidence for structure in the ejecta from the asymmetry in the [Ne II] 
and [O IV] lines. 

Thus, a lower limit to the abundance of neon in QU Vul with respect to the Solar value 
given by the strength of a given neon forbidden line alone is: 



L H J QUVul ^ 

—rwri - 



Nh 



7.413 X 10-5 



(6) 



where the Solar abundance of neon by number [Ne/H] = 7.413 x 10-^ flLoddersll2003l ). and 



53 



the total number of hydrogen atoms in the nova shell is taken as {M^j/nip) = 4.3 x 10 
On dates when more than one neon forbidden line was observed, the individual abundances 
derived from eqn. [6] are summed to provide an estimate of the total neon abundance. We 



^http://www. physics. nist.gov 
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stress that the total mass of neon in the ejecta of QU Vul, derived by the method described 
above, is a lower limit. The oxygen abundance is calc ulated in a sim ilar fashion, where the 



solar abundance of oxygen is [O / H]^ = 4.898 x 10 (lLoddersll2003l ) 



The observed neon line flux densities, log(Te 



level populations, and derived neon 



abundances are summarized in Table [T] An initial value of log(Te) = 4.0 was adopted for 
the abundance calculation for the first few hundred days. Later, a value of log(Te) = 3.0 was 
adopted based on the continuum-subtracted hydrogen Brackett-7 to Paschen-a emission line 
ratio derived from the Hubbl e Spac e Telescope (HST) narrow-band observations of the ejecta 
obtained bv iKrautter et al.l J2002h on day 2783.9 (= JD 2,448,843.5). Assuming Case-B 
condit ions and the hydrogen recombination line emissivities tabulated by lStorey &: Hummer 
f ll995h . the HST [HI{7 - A)/HI{A - 3)] = 7.27 ± 0.46 x lO^^ ^atio suggests that log(re) = 



3.0 for a ng = 1.0 x 10^ cm ^. The latter rig is of the order expected from eqnO 

Examination of the line profiles in the optical spectrum (Fig. |2]) that we obtained ~ 
242 days prior to the Spitzer observations lead us to conclude that the Balmer, He I, and 
He II lines probably come from both the low temperature ejecta and from the hot accretion 
disk in the binary system. Given the resolution of our spectrum, it is impossible to de- 
convolve the different components, and an analysis of the hydrogen line ratios yields an 
upper limit to the electron density in the ejecta. Dereddened (assuming E(B-V) = 0.6, 
Andrea et al.lll994l : ISaizar et al.lll992l ) flux ratios derived from Gaussian fits to the prominent 
H recombination lines in the day 6836 (= JD 2,452,895.8) spectrum yield {Hf3/Ha) ~ 
0.29 ± 0.05 and {Hj/Ha) ~ 0.14 ± 0.03. Assuming Case-B c ondit ions and the hydrogen 
recombination line emissivities tabulated by lStorey fc Hummer! (119951 ). the optical line ratios 
are consistent with log(Te) = 3.0 and He ~ 200 — 1,000 cm~^. Again, this range of values 
for Tie is consistent with the the value we derive by applying eqn [51 

For comparison, we have re-calc ulated the neon ab undances that can be i nferred from 



the d ata taken on 1985 May 15.4 UT flGehrz et al.lll985r). 1985 August 23.4 U T flGehrz et al. 



19861 ). and 1986 July 30.0 UT and August 11.0 UT flGreenhouse et al.lll988h based upon a 
more recent understanding of the expansion velocity and density development of the ejecta 
as described above, as well as the new values for the excitation and emission constants 
associated with the observed neon transitions. As is evident from the data in Table [H the 
early epochs clearly fail the criterion that Ue << ficrit, suggesting that the dramatic difference 
in abundance estimates for the early epoch relative to the later epochs likely is due to strong 
collisional damping of the emission lines. Also, it is not obvious that the relative locations 
of line formation within the expanding ejecta have remained the same over time. 

The large neon overabundances (at least 76 times solar on day 2783.9 and 168 times 
solar on day 7077.6) derived from our most recent observations (Table 1, column 13) are 
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even larger than the value of 21.7 ± 1.7 reported by Schwarz (2002; hereafter S02) based on 
detailed photoionization models constrained by multi-wavelength line fluxes obtained much 
earlier in the development of the ejecta. Closer examination of the assumptions made in the 
current analysis reveals that the discrepancies between the current analysis and that made 
by S02 can be easily resolved. The best fit CLOUDY model parameters reported by S02 
suggest a distance estimate of 2.4 kpc compared to our adopted value of 3.14 kpc. From 
equations (1), (4), and the the relationship between flux and luminosity, the abundance 
estimates scale as D"^. If we adopt instead a distance of 2.4 kpc, our calculated abundances 
drop by a factor of ~ 1.7. Further, in our analysis, we have made the elementary assumption 
that the ejecta uniformly fills a spherical volume with R = Vq At. If instead we assume 
a spherically symmetric shell with a total volume half that of the sphere approximation, 
the electron number density would go up by a factor of two consequently reducing our 
abundance estimates by a corresponding factor of two, bringing them into close agreement 
with estimates made by S02. On the other hand, we note that S02 used a solar neon 
abundance 1.59 higher than the more current value given by Lodders (2003). Correcting 
S02's values upward by this factor brings them in line with our results given the distance 
and the volume estimates discussed above. We also stress that the abundance values for neon 
given by us are inferred from the total mass of hydrogen in the ejecta as determined from 
radio continuum measurements (Taylor et al. 1988) and not from hydrogen recombination 
line strengths. 

Other estimates of the neon abundance by number in QU Vul using ionization correction 
factor (ICF) techniques and CLOUDY modeling have ranged from 60 - 254 times solar. 
The weaknesses in the ICF method and other CLOUDY models are discussed in detail 
by S02. The reader should be aware that Livio & Truran (1994) have pointed out that 
the uncertainties in the various methods of determining abundances in novae can lead to 
discrepancies of factors of ~ 2-4 under some circumstances. Our results (Table 1), taken in 
context along side rigorous photoionization modelling, suggest a significant overabundance 
of neon in the ejecta of QU Vul. 

The neon that we have detected in our spectra is most likely ^°Ne dredged-up from the 
WD outer layers at some time during the thermonuclear runaway. In addition, a fraction 
could be ^^Ne which is produced from the decay of ^^Na {ti/2 = 2.6 yr) and which is pre - 



dicted to be produced by thermonuclear runaways on ONeMg WDs (jStarrfield et al.lll992l ). 



Moreover, the formation and ejection of low levels of ^^Na, decaying into ^^Ne, can a ccount 



for th e production of some of the ^^Ne {Neon — E) anomalies found in meteorites (IBlack 



1972|). 
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4. CONCLUSIONS 



Given the data summarized in Tabled], we conclude that neon in this nova was at least 76 
and 168 times overabundant with respect to the solar value at epochs of 7.62 and 19.38 years 
past outburst. While this indicates that individual novae such as QU Vul might be important 
local sources for enriching the interstellar medium with neon, they likely contribute no more 
than a few percent of the total Galactic abundance of neon. However, a fraction of the 
ejected neon could be ^^Ne produced by the decay of ^^Na which may be produced during 
the outburst and could contribute to the ^^Ne {Neon — E) anomalies found in meteorites. 
However, the most significant implication of the high level of neon enrichment observed in QU 
Vul - since its source must certainly be dredge-up of matter from the underlying degenerate 
core - is that it provides unambiguous evidence for the occurrence of ONeMg white dwarfs 
among cataclysmic variable systems. The level of enrichment seen here is consistent with 
the high values pre dicted by the theoretical models of TNRs on ONeMg WDs presented by 



Ihadis et all tOQ'I }. 



We thank J. Jose and K. Vanlandingham for illuminating discussions. This work is 
based in part on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under NASA contract 
1407. RDG, CEW, LAH, and EFP are supported in part by NASA through contracts 
1256406, 1215746, and 1267992 issued by JPL/Caltech to the University of Minnesota. SS 
acknowledges partial support to ASU from the NSF and NASA. JWT is supported in part 
by the NSF under Grant PHY 002-16783 for the Physics Frontier Center, Joint Institute for 
Nuclear Astrophysics and by the DOE, Office of Nuclear Physics under Grant W-31-1009- 
ENG-38. 

Facilities: 

Spitzer (IRS), 

200 inch Hale Telescope (SpectroCam-10) 
MMT (Blue channel spectrograph) 



- 11 - 



REFERENCES 

Andrea, J., Derschel, H. & Starrfield, S. 1994, A&A, 291, 869 

Bevington, P. R., & Robinson, D. K., 1992, Data Reduction and Error Analysis for the 
Physical Sciences, McGraw-Hill: New York 

Black, D. C. 1972, Geochim. Cosmochim. Acta, 36, 377 

Blum, R. D., & Pradhan, A. K. 1992, ApJS, 80, 425 

Cohen, M., Walker, R. G., Barlow, M. J., & Deacon, J. R. 1992, AJ, 104, 1650 

Ferland, G. J. 1998, in "Wild Stars In The Old West: Proceedings of the 13th North Amer- 
ican Workshop on Cataclysmic Variables and Related Objects", eds. S. Howell, E. 
Kuulkers, and C. Woodward, ASP, 137, 165 

Gehrz, R. D. 1988, ARA&A, 26, 377. 

Gehrz, R. D. 1999, in Processes in Astrophysical Fluids,, eds.: O. Regev and D. Prialnik, 
Physics Reports, 311, 405 

Gehrz, R. D. 2002, in Classical Nova Explosions, eds.: M. Hernanz and J. Jose, AIP Con- 
ference Proceedings, Vol. 637, AIP: Melville, New York, p. 198 

Gehrz, R. D. 2008, in Classical Novae, second edition, eds.: M. F. Bode & A. Evans, Cam- 
bridge University Press, Cambridge, Chapter 8 

Gehrz, R. D., Truran, J. W., Williams, R. E., & Starrfield, S. G. 1998, PASP, 110, 3 

Gehrz, R. D., Grasdalen, G. L., & Hackwell, J. A. 1985 ApJ, 298 L163: erratum 1986, ApJ, 
306, L49 

Gehrz, R. D., Grasdalen, G. L., Greenhouse, M. A., Hackwell, J. A., Hayward, T., & Bentley, 
A. F. 1986, ApJ, 308, L63 

Gehrz, R. D. et al. 2007, Rev. Sci. lustrum., 78, 011302. 

Greenhouse, M. A., Grasdalen, G. L., Hayward, T. L., Gehrz, R. D., & Jones, T. J. 1988, 
AJ, 95, 172 

Hayward, T. L., Miles, J. W., Houck, J. R., Gull, G. E., & Schoenwald, J. 1993, Proc. SPIE, 
1946, 334 

Houck, J., et al. 2004, ApJS, 154, 10 



- 12 - 



Hummer, D. G., Berrington, K. A., Eissner, W., Pradhan, A. K., Saraph, H. E., & TuUy, 
J. A. 1993,A&A, 279, 298 

Iliadis, C, Champagne, A., Jose, J., Starrfield, S., & Tupper, P. 2002, ApJS, 142, 105 

Jones, B., & Puetter, R. 1993, Proc. SPIE, 1946, 610 

Kramida, A. E, & Nave, G. 2006, Eur. Phys. J., 37, 1 

Krautter, J., Woodward, C. E., et al. 2002, AJ, 124, 2888 

Livio, M., & Truran, J. W. 1994, ApJ, 425, 797 

Lodders, K. 2003, ApJ, 591, 1220 

McLaughlin, B. M. & Bell, K. L. 2000, J. Phys. B., 33. 597 

Osterbrock, D. E., 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei, Mill 
Valley: University Science Books 

Pradhan, A. K. & Peng, J. 1995, in The Analysis of Emission Lines, eds.: Robert E. Williams 
and Mario Livio, Cambridge Univ. Press: New York, p. 8 

Rosino, G., et al. 2002, A&A, 257, 603 

Saraph, H. E. & TuUy, J. A. 1994, A&AS, 107, 29 

Saizar, P. et al. 1992, ApJ, 398, 651 

Schmidt, G. D., Weymann, R. J., & Foltz, C. B. 1989, PASP, 101, 713 

Schwarz, G. J. 2002, ApJ, 577, 940 (S02) 

Starrfield, S. et al. 1992, ApJ, 391, L71 

Storey P. J. & Hummer, D. G. 1995, MNRAS, 272, 41 

Taylor, A. R., Hjellming, R. M., Seaquist, E. R., & Gehrz, R. D. 1988, Nature, 335, 235 
Werner, M. W., et al. 2004, ApJS, 154, 1 



This preprint was prepared with the A AS macros v5.2. 



-13- 



Table 1. Physical Parameters for Neon Emission Lines Observed in QU Vulpeculae 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


_^(12)- 








Emission 


Integrated 














[tj]qUVuI 




Date 




Line 


Flux 


log Te 






"^crit 








OBS 


(UT) 


Day 


(Aim) 


(10-12 grg s-i cm-2) 


(K) 


T ■ ■ 


(10^ cm-3) 


(10^ cm-^) 


(10™ 


) (10-*9) 


[#]0 


WIRO 


1985 May 15.4 


141.3 


[Ne II] 12.81 


220±20 


4.0 


0.283 


332 


7.03 


4.5 


19.5 


^20 


WIRO 


1985 Aug 23.4 


241.3 


[Ne II] 12.81 


97±10 


4.0 


0.283 


67 


7.03 


2.1 


8.6 


^9.4 


IRTF 


1986 Jul 30.0 


581.9 


[Nc II] 12.81 


71.9±3.4 


4.0 


0.283 


4.8 


7.03 


3.5 


6.4 


^13 


IRTF 


1986 Jul 30.0 


581.9 


[Ne VI] 7.64 


247±19 


4.0 


2.720 


4.8 


3.44 


0.58 


5.5 


^3.6 


WIRO 


1986 Aug 11.0 


593.0 


[Ne VI] 7.64 


249±25 


4.0 


2.720 


4.8 


3.44 


0.58 


5.6 


^3.6 


HALE 


1992 Aug 09.3 


2783.9 


[Ne II] 12.81 


6.84±1.09 


3.0 


0.272 


0.04 


2.31 


24.1 


0.61 


^76 


SST 


2004 May 11.7 


7077.6 


[Nc II] 12.81 


0.475±0.032 ^ 


3.0 


0.272 


0.0026 


2.31 


23.1 


0.04 


^73 


SST 


2004 May 11.7 


7077.6 


[Nc 111] 15.56 


1.09±0.028 


3.0 


0.596 


0.0026 


1.08 


30.1 


0.17 


^95 


SST 


2004 May 11.7 


7077.6 


Total Neon 
















^168 


SST 


2004 May 11.7 


7077.6 


[O IV] 25.91 


0.304±0.012 


3.0 


1.641 


0.0026 


0.02 


4.9 


0.001 


^2.3 



^Column 12 represents the abundance by number of the given emission line. 

''The IRS 12.81 /im values are the weighted mean of measurements from the SL and SH modules. 
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wavelength (//.m) 

Fig. 1.— SpcctroCam-10 spectrum of QU Vul taken on 1992 August 09.3 UT. The [Ne II] 
line at 12.81 iim has an integrated flux of 6.84 x 10~^^ erg s~^ cm~^. 
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Fig. 2. — MMT Blue Channel optical spectrum of QU Vul taken on 2003 September 
13.31 UT. Prominent recombination lines of H and He are evident. The [Ne III] lines near 
3869 and 3968 A are weak and blended with adjacent hydrogen lines at the spectral resolution 
of the instrument. 
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Fig. 3. — Spitzer background-subtracted Short Low (SL) spectrum of QU Vul taken on 2004 
May 11.7 UT. The [Ne 11] line at 12.81 /im has an integrated line flux o£ 4.96 x 10"" erg 

— 1 —2 

s cm ^. 
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Fig. 4.— ^piteer Short High (SH) spectrum of QU Vul taken on 2004 May 11.7 UT. No 

sky subtraction was performed. The [Ne II] hne at 12.81 //m has an integrated line flux of 
4.34 X 10"^^ erg s^^ cm^^. Also evident is the [Ne III] line at 15.56 /im with an integrated 
line flux of 1.09 x 10~^^ erg s~^ cm~^. For clarity, a representative subset of continuum errors 
is plotted. 
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Fig. 5. — Spitzer Long High 1 (LHl) spectrum of QU Vul taken on 2004 May 11.7 UT. No 
sky subtraction was performed. The [O IV] hne at 25.91 /xm has an integrated hne flux of 
3.04 X 10~^^ erg cm~^. The absolute flux cahbration of the flux near the line is uncertain 
by less than 10%. 



